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Abstract 

Orbital detectors without pointing capability have to keep their field of view 
axis laying on their orbital plane, to observe the largest sky fraction. A general 
approach to estimate the exposure of each sky element for such detectors is a valu- 
able tool in the R&D phase of a project, when the detector characteristics are still 
to be fixed. An analytical method to estimate the sky exposure is developed, which 
makes only few very reasonable approximations. The formulae obtained with this 
method are used to compute the histogram of the sky exposure of a hypothetical 
gamma-ray detector installed on the ISS. The C-n- code used in this example is 
freely available on the http : // cern . ch/ casadei/ softw are . html web page. 

1 Introduction 

In addition to the large, complex and very expensive orbiting observatories, smaller size 
and lower cost detectors, whose goals are restricted to a narrower field, can perform in- 
teresting astrophysical measurements. To be (relatively) cheap, a satellite cannot have 
special pointing capabilities, beyond the normal gyroscopic stabilization and control, 
which is necessary not to loose communication with Earth stations and to keep safe 
satellite orientation. 

When trying to define what are the parameters which should characterize an or- 
biting detector, one would like to be able to get values which are good estimations of 
the final parameters much before developing a complete Monte Carlo simulation. As 
an example, let us consider a detector which aims to measure the cosmic gamma-ray 
background in a wide energy range: its design acceptance and the mission duration 
will depend on the expected flux of known sources, in addition to the detector energy 
resolution. 

To estimate the detector field of view and sensitive area, together with the mission 
duration, one needs to know the measured flux (per energy bin) coming from each sky 
element, which depends on the sources and on the detector characteristics. Thus, one 
will first produce the histograms of the sky exposure for different mission durations 
and then will consider different sensitive areas (i.e. different detector configurations), 
in order to find the configuration which best matches the expected fluxes. 

Another case in which a method to determine the sky exposure for a given detecor 
can be very useful is when one tries to understand the possible performance of an 
existing detector for the search for a new "channel" in the data analysis, and the access 
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to the detailed simulation software is not easy, fast, or allowed. If the obtained values 
suggest the possibility of doing interesting astrophysics, one will of course use the 
official detector simulation program to refine the estimation. 

In this paper, an analytical method to obtain the sky coverage of an orbiting detector 
is shown, which makes use of few reasonable assumptions to simplify the computation: 
(1) the orbit eccentricity is zero; (2) the orbit precession period is much larger than 
the revolution period; (3) the detector field of view is a cone, whose axis lies in the 
orbital plane, rotating with the same angular velocity of the satellite revolution. These 
assumptions are valid in the case of a detector installed on the International Space 
Station' (ISS), taken as an example in the following sections. A program has been 
developed which makes use of this method: its C++ code is freely available on the 
ihttp : / / cern . ch/ casadei/sof tware . html web page. 

2 The method assumptions and algorithm 

2.1 Hypotheses 

The first hypothesis (null eccentricity e) implies that the revolution velocity is constant. 
Its explicit dependence on the orbit parameters is shown in ref. yj: the result is a 
function with a rather weak dependence on e, hence this approximation is quite good 
for small values of e. 

The second hypothesis (orbit precession period T » t, where t is the revolution 
period) is valid in general, and allows for a useful approximation: the method first 
computes the sky exposure over a single orbit, then rotates this map following the orbit 
precession. 

The hypothesis that the field of view is a cone with axis in the orbital plane is not 
true in general. The real detector acceptance is usally different from a uniform cone, 
though one often accepts tracks within some "fiducial" cone to make easier the off-line 
data analysis (avoiding "border effects"). Having the axis contained in the orbital plane 
maximizes the fraction of the sky seen by the detector, hence it is a quite natural choice. 
Finally, having angular velocity equal to that of the revolution motion means that the 
angle (in the orbital plane) between the cone axis and the radial vector (from the orbit 
center to the satellite) is kept constant. This is true for detectors installed on the ISS. 

The detector acceptance will in general depend both on the incident angle a of the 
photon and on its energy, with a - Q for tracks parallel to the cone axis. However, 
we neglect this angular dependency and consider a uniform acceptance in a inside a 
cone with half-aperture a^. The analytical method presented here is really fast^ from 
the computing point of view, and it is easy to change the source code to implement the 
effect of the a dependence of the acceptance (only a single file needs to be changed). 
On the other hand, it will be much more tricky to implement the azimuthal dependence 
of the detector acceptance: if the user really needs to reach this level of detail, she/he 
would benefit from switching to the official detector simulation. 

' http : / /www . hq . nasa . gov/os£/station/viewing/issvis . html 

^To compute the final histogram of this paper, shown in figure^ a Linux based PC with a Pentium III 1 
GHz CPU takes 1.8 s. 



2 



2.2 Reference systems 



In this paper, three reference systems are used: the orbit reference system (hereafter 
ORS), the galactic reference system (GRS), and the precession reference system (PRS). 
The GRS is the system in which we usually want to express the results. The PRS is 
used as intermediate step by the algorithm, to simplify the precession computation. 
Finally, the equatorial (x, y) plane of the ORS is the plane of the orbit. PRS and GRS 
are inertial reference sytems (they do not rotate). 

Because we will consider here only mission durations which are at least few times 
longer than the precession period-^, the exact position of the detector when it starts 
data taking is not important. Hence, we are only interested in the orbit unit vector 
h, which defines the z axis in the ORS. The {x,y) plane is the orbital plane in this 
reference system, and we are free to choose the axes direction in ordert to simplify the 
computation. 

A possible choice for the reference axes {X, Y, Z) in the GRS is the following. The 
axes are defined by the three unit vectors /, / and K respectively: / points to the vernal 
equinox direction, K to the North Pole and / completes the right-handed triplet | H. 
Any point in the sky can be addressed in the GRS using two angles: e [0,7r], which 
is the angle formed with K, and O e [0, 27r], expressing the rotation in the equatorial 
(X, Y) plane, starting from the / unit vector Astronomers usually adopt the "equatorial 
coordinates" declination S - n/l - @ and right ascension a - <S> instead of and 
However, using the colatitude instead of 6 simplifies the histogramming. 

The orbit precession is described by the constant unit vector k about which the orbit 
normal h rotates while keeping k-h(t) - cos /3 constant. The k unit vector defines the z' 
axis of the PRS, and we are free to choose the direction of the x' and y' axes to simplify 
the computation. The algorithm needs the PRS when computing the orbit precession 
eff'ect on the sky exposure map, which in the PRS is simply a rotation around k of the 
histogram obtained for a single orbit. 

2.3 Algorithm 

For each sky element, its exposure is equal to the product of its solid angle and the total 
time for which it is found inside the detector acceptance. During few orbits, as long as 
dh/ dt can be neglected, the same points are seen by the detector few times: if < 90 
deg (as usual), the celestial sphere is divided in two parts by the "belt" viewed by the 
detector. The exposure is largest for points situated along the central part of the belt, 
and decreases to zero going from the center to the boundaries. 

After a large number of orbits, the precession of the belt increases the sky fraction 
viewed by the detector. For each sky element, its exposure is equal to the convolution 
between the time spent inside the belt (which rotates following the orbit precession) 
and the exposure function inside the belt (computed along one revolution). 

The algorithm, for the case in which the detector acceptance is uniform inside a 
cone with half-aperture whose axis is defined as the detector axis, is based on the 
following ideas: 

1 . to find the exposure map of the sky in the ORS for one revolution, we will follow 
an analytical treatment (section|3|i; 

^Increasing somewhat the mission duration is simpler and cheaper than increasing the detector sensitive 
area (i.e. the pay load weight). 
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Figure 1: Sky sphere viewed by a cone. 



2. the (slow) orbit precession in the PRS is computed with a discrete algorithm 
acting on the exposure histogram (section lOl : 

3. the exposure map of the sky in the GRS is found analytically ("section l4!2t . 

If the real detector acceptance does not have cylindrical symmetry, one has first to 
compute the real belt exposure along one orbit by changing step 1 (i.e. by changing 
a single C function in the code), and then can proceed as before. When considering 
the dependence of the acceptance on the photon energy, one will repeat the three steps 
independently in each energy bin. 

3 Single orbit exposure map 

Let the detector field of view be a cone of half-aperture (figure O, with axis along 
the unit vector r (the "local vertical")'*. The latter turns around h with angular velocity 
o) - 2n/T, and instant velocity unit vector v = hxr (the "local horizontal")^. 

3.1 Sky coverage in the orbit reference system 

All points along the r trajectory are inside the cone for a time f = 2ac/cL) = q-ct/tt during 
each orbit. The fraction of time during which a point P in the sky is found inside the 
cone depends on the angle a formed with the equatorial plane (which contains r and 

^If the detector axis, which lies in the (x, y) plane, forms a constant angle with f, the result is exactly the 
same. 

'The three unit vectors n, v and r, form the basis of the so-called "local-vertical local-horizontal" (LVLH) 
reference system. 
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Figure 2: Side view of the sky sphere. 



Figure 3: Top view of the sky sphere. 



v). This time t(a) is maximum when a — and goes to zero when a approaches a^. 
For higher angles, the point is never contained by the cone (we restrict our analysis to 
ffc < 90 deg). 

The "parallel" passing through P is a circumference with length i = IttR cos a 
(figure IS. Let t// be the half angular width of this parallel, when projected to the 
equatorial plane (figure^: when a approaches ac, 4' goes to zero; as a goes to zero, ij/ 
approaches a^- The exposure time t{a) is equal to the ratio between the parallel width 
lij/ and the angular velocity: 

2 Ma) T 

f(a)= = -iA(a). (1) 

(Jj 71 

In order to find (//(a), we consider the triangle ABO on the equatorial plane (see 
figure and observe that the segment BO, the radius of the parallel passing through 
P, is r = /? cos a. The lengths of the segments AO and AB are R cos ac and r sin i/r 
respectively. The Pitagora's theorem implies that 

cos^ a - cos^ Cc + cos^ a sin^ i// , (2) 

hence one finds that sin if/ can be written as: 



sintA = — . (3) 

V cos^ a 

One can easily verify that this expression has the expected behavior when a approaches 
zero or a^. 

Finally, the time t(a) during which the point P, seen at the angle a with respect to 
the equatorial plane, is viewed by the detector is: 



t(a) - - arcsin -v/ 1 - " ., , (4) 
n y cos^ a 

with -Oc < a < (see figure^] where = 25 deg). 

A better choice of the angle is the orbital colatitude = | - a (also shown in 
figure^}: for any direction d in the sky one has h-d - cos 0. Every sky element can be 
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Figure 4: Fractional exposure time (in units of the orbit period t) of a cone with 
ffic - 25 deg, as function of the angles a and 6. 



mapped using and the angle (f> of rotation^ about h: the solid angle element centered 
on d is dO. = sin 9 d6 d(p and its exposure per orbit is t{d) dQ., i.e. : 



^ T I COS^ Cf, 

t(d) dQ.- - arcsin -v / 1 - 
n V 



sin^e 



sin6ld6'd(^ 



(5) 



for (|-Q'c)<0<(|+Q'c) and f(d') dQ = otherwise. In the following, we will omit 
to recall that the exposure is zero when the condition (| - a^) < < (j + a^) is not 
satisfied. 

In numeric computations, given the partition of the sky in elements AO. = A(cos 0) A<p, 
corresponding to uniform binning in (p and cos 0, one can produce the exposure his- 
togram using equation (|5jl, which can be rewritten as: 



^ T I COS (Y 

t{d) AQ = - arcsin -i/ 1 ^ A0 A(cos 0) , 

n V 1 - cos^ 

where cos is the central value of the considered bin, whose width is A(cos 0). 



(6) 



3.2 Sky coverage in the precession reference system 

Equation (|6j is valid in any (x, y, z) reference system where h is the unit vector defining 
the z-axis. In order to compute the effect of the orbit precession, we move from the 
ORS to the (jt:',y,z') PRS with a rotation about the x - x! axis^; 



A = 



(10 
cosy6 sin/3 
-sinyS cosy6^ 



(7) 



In the PRS the orbit precession simply results in a drift along 0' (<p' - along the x' 
direction, increasing following the orbital motion) of the sky elements in the cos 6*' vs. 

= define the direction of the x axis. 
'We are free to choose the direction of both x and x' axes, because the result is independent from this 
choice. 
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Figure 5: Single orbit exposure map per unit period, for a detector with field of view 
of 50 deg (25 deg is the cone half aperture), orbiting on board of the the ISS. The 
ExposureFraction histogram size is 100 x 100 bins. 

<p' exposure histogram (i.e. in a horizontal shift of the column contents). We consider 
the sky element A<p'A{cos 9') centered around the direction which can be written in 
Cartesian coordinates as 



z' - cos 6' 

To apply equation (|6j, which is independent of (f>, we only need to find z - cos by 
applying A obtaining: 



Using equation together with (|6}, one can obtain a histogram similar to figure|5] 
which shows the sky exposure map (for one orbit, in units of the revolution period) seen 
by a detector whose field of view is a cone with half aperture 25 deg, installed on the 
ISS: in the PRS, the "belt" is inclined of 51.6 deg, which is the ISS orbit inchnation^. 

*The k unit vector points to the North celestial pole, hence one may use the celestial coordinates (decli- 
nation and right ascension) in the PRS. This choice fixes the direction of the x = x' axis. 



x' - sin 6' cos 4>' 
y' - sin ff sin (p' 



(8) 



cos 6 - z- Vl - cos^ 0' sin 4>' sin p + cos 6' cos (3 . 



(9) 
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In the following, we will illustrate the method by choosing a mission duration T^i - 
1000 days and the detector axis to be parallel to the ISS local vertical. These parameters 
are similar (though not identical) to those of the AMS-02 detector |2], which will be 
installed on the ISS in the near future. The ISS orbit parameters are approximated as 
follows: the revolution period r ^ 1.5 h, the precession period Tpr ^ 2 months. With 
altitudes in the range 350^00 km above sea level, the orbit eccentricity e - 6.55 x lO"'* 
is very small (i.e. our first hypothesis is good). In addition, the ISS attitude is controlled 
in order to keep its local vertical axis parallel to the radial vector (hence the third 
hypothesis is also good). 

4 Orbit precession 

Let us call ExposureFraction the histogram representing the "belt" in the PRS (fig- 
ure|5}, and PrsSkyMap the histogram of the exposure in the PRS, which we have now 
to fill. The bin contents of ExposureFraction are pure numbers: they must be mul- 
tiplied by the orbit period r to obtain the exposure of each Af2 sky element during a 
single orbit. 

4.1 Orbit precession in the precession reference system 

To fill the PrsSkyMap histogram, which has A^^ bins in 0' and Ne bins in cos0', we 
have to compute the effects of the orbit precession, which in the PRS results in a drift 
along (p' of the sky elements obtained using equation (|6j after the substitution The 
algorithm is the following. 

1 . One has first to compute the time required for the orbit precession to move each 
bin on top of the following one. This time is - TpJN^, with Tpi = 2 months is 
the precession period. 

2. Given the mission duration T^i, one has to find the number of times the ^ bins 
have to be rotated, which is A^i-o - T^i/t^. 

3. Then one has to performs A^io iterations which consist of a loop over the A^^ 
columns, whose single step is the sum of the ExposureFraction (/- l)-th column 
contents to the PrsSkylVlap i-th column (with / decreasing from A^^ to 1). Of 
course, the last column should be added to the first one. 

4. Finally, to get the exposure time for each AQ sky element, one has to multiply 
each PrsSkylVlap bin by the revolution period t. 

Figure|6lshows the result for our hypothetical detector on the ISS. The PrsSkyMap 
histogram has 100 x 100 bins. 

4.2 Orbit precession in the galactic reference system 

The histogram GrsSkyMap representing the sky map in the GRS can be filled with 
the assumption that the O-granularity of GrsSkyMap should be identical to the (f>'- 
granularity of PrsSkyMap (which is quite reasonable). However, one can not find the 
PRS histogram and then transform it to the GRS map, because continuous transforma- 
tions can not be applied to discretized maps without introducing distortions. Then the 
procedure is the following. 
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Figure 6: Exposure (in minutes) in the PRS of a detector with field of view of 50 deg 
(25 deg is the cone half aperture), orbiting on board of the the ISS. The orbit period is 
90 min, the precession period is 2 months, and the mission duration is 1000 days. The 
PrsSkyMap histogram size is 100 x 100 bins. 



For each GrsSkyMap bin, relative to the direction d corresponding to the coordi- 
nates (3), cos 0) of the bin center, we need to find the PRS coordinates which represent 
the same direction: 

dms = (0'(0, <!•), cos 6i'(0, 0)) = (0}, cos 6K) , (10) 

where i g [0,Ng - 1] and j g [0,A'^ - 1] are the row and column indices, respectively, 
of the PrsSkyMap bin which contains rfpRs- 

Each ^'-bin (and O-bin, for the previous assumption) is A<p' = In/N^ wide and is 
inside the detector field of view for a time per orbit. In the PRS, the time T{d) spent 
inside this bin is the sum over Nro rotations of the PrsSkyMap columns, due to the orbit 
precession. Hence 

r(d)= Varcsin, 1--^^^^ (11) 
V 1- W,0') 

with 

Fkiff, <p') = cos^ e{&i, <p'j - kA(f>') = cos^ 0(0;, 0;.(^)) (12) 
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Figure 7: Exposure (in minutes) in the GRS, after an Euler's rotation of angles (-33, 
62.6, -77.75) deg. 

where jik) is the ^'-bm which contains (f)'j - kh(p' . Finally, to obtain the total expo- 
sure in the considered GrsSkyMap bin, we have to multiply the time spent viewing 
each sky portion centered around (<1>, cos ©) for A^io times its solid angle: the result is 
Ni-oT{d)AQ., where T{d) is given by the formula il 1> . 

Actually, the real algorithm follow by the program starts from the GrsSkyMap 
histogram and loops over its bins. For each bin in the GRS, it finds the corresponding 
bin in the PRS using an Euler's rotation**, then applies formula il 1\ . Figure0shows the 
final result, after a rotation with Euler's angles of (-33, 62.6, -77.75) deg. The first angle 
is the rotation about the third axis, which transform the first axis in the 'line of nodes'. 
The second angle is a rotation about the line of nodes, which moves the third axis in the 
final position. The third angle is a rotation about the new third axis, which transforms 
the line of nodes into the final first axis. In order to show another possible choice of the 
coordinate system, we considered the PRS as the celestial-equatorial reference system 
in our example, and chose the Euler's angles which transform it to the GRS, expressed 
in the B1950.0 epoch as defined in 1958 by the lAU H. 

;http : //mathwor Id . wolfram . com/Euler Angles . html| 
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5 Detector sensitivity 



Having obtained the exposure time of each sky element, the detector sensitivity is found 
in the following way. The first step is to compute the product of the effective detector 
sensitive area A(E) (usually expressed in cm^) with the sky element exposure AQ f (sr 
s): 

S (E) = A{E) AQ t (cm^ sr s) . (13) 

The effective area A(E) in general will depend on the photon energy E and on the 
incident direction. Here we consider a uniform acceptance over all directions inside 
the field of view, hence our expression for S (E) is equivalent to the average over all 
directions of the real detector efficiency. 

If a known source, covering the sky element AQ (which can be contained in a 
single histogram bin or extend over a set of adjacent bins), has a photon flux F(E) 
(cm^^sr"' s"') in a particular energy bin centered at E, then the number of events 
measured by the detector are expected to be: 

N,^e.AE) = F{E)S(E). (14) 

On the other hand, there will be some background flux of photons coming from the 
same sky element, which will produce background events in each energy bin. Formula 
il4\ can also be used to find the number of expected background hits in the detector. In 
case one considers a single energy bin, she/he would define the detector source sensi- 
tivity in this bin as the ratio between the number of expected measured events and the 
square root of the number of background hits'", or follow the more detailed approach 
of ref. f^. However, the usual approach is to consider the whole detector energy range 
at the same time, and to perform a fit of the measured distribution (deconvolved with 
the detector sensitivity) with a function which is modeling both the source spectrum 
and the background flux at the same time. 

In this way, for each set of detector parameters, it is possible to identify the sources 
which are visible with the considered orbit and mission duration. 

6 Conclusion 

The analytical method developed in this paper allows for a fast estimation of the detec- 
tor sensitivity to known sources, when the hypotheses considered in section|2are valid. 
The source code, freely available on the author's web site", is a useful tool during the 
design phase of an experiment making use of an orbital detector, when the goal is the 
definition of its parameters and of the mission duration. In addition, it is useful also 
as a way to easily check the performance of existing detectors, when access to their 
detailed simulation program is not easy or allowed. 

On the other hand, if one really needs a precise estimation of the detector perfor- 
mance, she/he would need to make use (or to develop) a program which exactly tracks 
the satellite along its orbital motion usingcomphcated formule, which can be found in 
papers like ref. 1 1 ] or textbooks like ref. \5] and f6] . In addition, the program should be 
able to correctly simulate the detector response to photons with different energies and 
incident directions, as was done for example by the AMS collaboration 
But the price for the higher precision is a CPU intensive and complex program (with 
usually harder user interface than this analytic method), and a long development time. 

'"The number of background hits per energy bin is assumed to be Poisson distributed. 
" http : //cern . ch/casadei/sof twar e . html| 
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